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In this paper we study Chaplygin's Reducibility Theorem and extend its applicability to nonholonomic 
systems with symmetry described by the Hamilton-Poincarc-d Alcmbert equations in arbitrary degrees of 
freedom. As special cases we extract the extension of the Theorem to nonholonomic Chaplygin systems 
with nonabelian symmetry groups as well as Euler-Poincare-Suslov systems in arbitrary degrees of free- 
dom. In the latter case, we also extend the Hamiltonization Theorem to nonholonomic systems which do 
not possess an invariant measure. Lastly, we extend previous work on conditionally variational systems 
using the results above. We illustrate the results through various examples of well-known nonholonomic 
systems. 



Introduction 

oo 

Although it is well known that nonholonomic mechanical systems are not variational [2] 
and thus their mechanics cannot be expressed in terms of canonical Hamilton equations, 
nevertheless several authors (dating back at least as early as S.A. Chaplygin [HI US] 
and Appell pQ) have attempted to express the mechanics of nonholonomic systems in 
Hamilton-like forms through several methods. Perhaps the most well-known of these 
methods is Chaplygin's own Reducibility Theorem, whose first part states that for non- 
holonomic systems in two generalized coordinates (qi, #2) possessing an invariant measure 
with density N(qi, 92), the equations of motion can be written in Hamiltonian form after 
the time reparameterization dr = Ndt (in this context iV is known as the reducing mul- 
tiplier, or simply the multiplier). The second part of the Theorem (see [19]) says that 
if a nonholonomic system can be written in Hamiltonian form after the time reparame- 
terization dr = f(q)dt, then the original system has an invariant measure with density 
f m ~ 1 {q), where m is the degrees of freedom and the function / is again called the reduc- 
ing multiplier, or simply the multiplier. Since both components of the theorem involve a 
reparameterization of a nonholonomic system into a Hamiltonian one, one often refers to 
this as the Hamiltonization of a nonholonomic system, although we shall refer to it here 
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as the Chaplygin Hamiltonization insteacf] 

Chaplygin's original motivation for such a Hamiltonization of nonholonomic systems seems 
to have been rooted in his interest in the explicit integrability of nonholonomic systems. 
Indeed, in [15] Chaplygin applies his method to integrate what would later become known 



as the Chaplygin sleigh (see Section 4.5), and remarks that his general procedure (using 
the reducing multiplier) for integrating certain two degree of freedom nonholonomic sys- 
tems is "interesting from a theoretical standpoint as a direct extension of the Jacobi 
method to simple nonholonomic systems." Chaplygin further applied his theorem to inte- 
grate other nonholonomic systems by quadrature [IE], as did Kharlamova later [29]. Thus, 
the reducing multiplier method has historically been interesting and important from the 
standpoint of the integrability of nonholonomic systems (for more historical notes on the 
origin of the theorem, see [37], or [18] for a more geometric viewpoint, and [8]). This mo- 
tivation led us to investigate extensions of the theorem in this work in the hope of further 
expanding its applications to the integrability of nonholonomic systems. One would like 
a more straightforward procedure which eliminates such guesswork. 

After the introduction of Chaplygin's theorem, subsequent research on the theorem has 
resulted in, among other things, an extension to the quasicoordinate context |34J, a study 
of the geometry behind the theorem [T8J EU 125] , discoveries of isomorphisms between 
nonholonomic systems through the use of the theorem [7J, an example of a system in 
higher dimensions Hamiltonizable through a similar time reparameterization [21], an in- 
vestigation of the necessary conditions for Hamiltonization for abelian Chaplygin systems 
[2"7] (see Section 1.1 for a definition), Poisson structures for rolling bodies without slip- 
ping [SHE]) and an investigation of rank two Poisson structures in nonholonomic systems 
[35] . In addition, the survey paper [8] presents, among other things, many of the known 
examples to which Chaplygin's theorem is applicable. However, two important aspects 
yet to be resolved are the extension of the theorem to general nonholonomic systems with 
symmetry of arbitrary degrees of freedom and, since the theorem rests on the availability 
of an invariant measure, we are also interested in applying a time reparameterization to 
"Hamiltonize" a nonholonomic system not possessing an invariant measure, where dif- 
ferent dynamical effects may arise P, El HD]- We note that in this case our resulting 
"Hamiltonization" of a system not possessing an invariant measure should perhaps more 
properly be called a "Poissonization," since it will in general result in a degenerate Poisson 
bracket satisfying the Jacobi identity (see Section 4.5 below), in addition to the continued 
non-existence of the invariant measure (however, we will continue to refer to this process 
as "Hamiltonization," keeping in mind this discussion). Moreover, the theorem is com- 
monly used in a rather guess-and-check manner, where one considers systems with known 
invariant measures and then guesses at the reducing multiplier based on the degrees of 
freedom. 

In this paper we consider the aforementioned questions for a general nonholonomic system 
with symmetry governed by the Hamilton-Poincare-d'Alembert equations. In Section 1 
we briefly discuss the mechanics of these systems, as well as for two special cases of them 
(nonabelian Chaplygin systems and Euler-Poincare-Suslov systems), and present results 
in Sections 2.1, 2.2 and 2.3 which generalize the first part of Chaplygin's Theorem to 



1 We introduce this term because there are other ways of writing the reduced, constrained mechanics of a 
nonholonomic system as a Hamiltonian system that do not involve a time reparameterization, for example as was 
done in I4ll23l. 



2 



these higher dimensional systems with symmetry, deriving the necessary conditions (in- 
dependent of the existence of an invariant measure) for Chaplygin Hamiltonization as a 
coupled set of first-order partial differential equations in /. These equations eliminate 
the guesswork discussed above, and in the special case of Euler-Poincare-Suslov systems, 
we present results in Section 2.3 which allow Chaplygin Hamiltonization even when the 
system does not posses an invariant measure. In Section 3 we use the previous results 
to extend the idea of conditionally variational systems introduced in [23] and apply it to 
Chaplygin Hamiltonize the entire nonholonomic system (reduced constrained equations 
plus the nonholonomic constraints). Lastly, we devote Section 4 to illustrating these 
results and showing how special cases of the results presented lead to some of the re- 
sults found in the works cited above, and discuss some relevant future directions in the 
Conclusion. 

1 Nonholonomic Systems with Symmetry 

Consider a nonholonomic system with an n dimensional configuration manifold Q and 
mechanical Lagrangian L which is subject to k linear nonholonomic constraints described 
by the distribution T> (moreover, we shall restrict our attention to mechanical Lagrangians 
for the remainder of the paper). Suppose that we have a Lie group G which acts freely 
and properly on the configuration space Q, with the Lagrangian L and constraints T> 
invariant with respect to the induced action of G on TQ. For simplicity, assume also that 
the constraints and the orbit directions span the entire tangent space to the configuration 
space: 

V q + T q Orb{q) = T q Q, 
sometimes known as the dimension assumption [2]. 

With sufficient regularity we can use the Legendre transform to pass to the constraint 
phase space M. = FL(V). The quotient space M. = M. /G is a smooth quotient manifold 
with projection map p : M. — > M., and all intrinsically defined vector fields then push down 
to A4, allowing one to write the equations of motion for the reduced constrained Hamil- 
tonian mechanics using a reduced almost-Poisson (in general) bracket on M.. The result- 
ing nonholonomic equations of motion are known as the Hamilton-Poincare-d'Alembert 
(HPD) equations and split into a coupled set of second-order equations on the shape 
space M := Q/G and first-order nonholonomic momentum equations on g* [2} [31], whose 
number equals s := dim S q , where S q := V q fl T g Orb(g). 

Following |2j we will now give some of the details, however before doing so let us fix the 
following index conventions. The indices a, b, c, . . . will range from 1 to k :— dim(Q) and 
correspond to the symmetry directions, i, j, . . . will range from 1 to s (s < k is the number 
of momentum equations) and correspond to the symmetry directions along the constraint 
space, and a,j3, . . . will represent the indices for the shape variable r e M := Q/G and 
range from 1 to m := n — k = dim(M) (the dimension of the shape space). Also, here and 
for the remainder the of the paper we shall enforce the Einstein summation convention, 
unless otherwise indicated. 
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Begin by constructing a body fixed basis eb{g,r) = Ad g e {r) as in [2], where g G G and 
r G M, such that the infinitesimal generators (ei(g,r))Q of its first s elements at a point q 
span iSq. Assuming G is a matrix group and ef is the component of ej(r) with respect to 
a fixed basis {b a } of the Lie algebra g, we can then represent the constraint distribution 
V as 

V = span{<$ef<V, -^fy. + d r «}, 

where we will denote by Q % the body angular velocity components of the constrained 
vertical space. 

Defining the induced coordinates (g a ,r a ,pi,p a ) on M. by 

Pi = 9dPae d i = Hdfl*, pa=Pa- fibA b a , 

where /i6g* and /i a are its components with respect to a fixed dual basis, p a = 8L/dg a 
and p a = dL/df a , the Hamilton-Poincare-d'Alembert equations on Ai are given by [21 [5]: 



~ _ ,, r a r bd 9h M dh M 

Pi - -/i a G M e i e i -^- + ^ a i'i j B-^-> ( l ) 

r " = % < 2 > 



along with the constraints 

Here hj^(r,Q,p) = pifl 1 + p a r a — Z c is the constrained reduced Hamiltonian (where 
l c (r,r,Q) = /(r, f,£ = — *4r + fie) is the constrained reduced Lagrangian), B% are the 
coefficients of the curvature of the nonholonomic connection: 

dA a dA a B 

™ a P_ _i_ /~ia AO Ac 

a P ~ dr? dr a hc ^ a ^ 

where C^ c are the structure constants of the Lie algebra g, £ b = (g~ 1 ) b a g a , and finally the 
Fife are given by 

-pa i I fia J> Ac 

W ~ ^bcZi-AfB- 

Moreover, in equations (0)- (J3]) the quantities /i a = dl/8C, a should be restricted to M. by 
substituting in the constraints Q. Hereafter we shall denote any expression into which 
the constraints have been substituted with a subscript c, as in (/x ) c . 
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Now, the equations ({T|-(|3| can be written with respect to an almost-Poisson (AP) bracket 

{•> -}m g iven b y 



dg dk dg dk 



dp a dp a dpi 

+ r 5 i— 

dr a dp f3 dppdr a ) P^Pf 3 dpadpp 



, dg dk dg dk\ dg dk 

+ i r >Pf>} +{P°>M7£T7£r> ( 5 ) 



where 



= -{^) c Cl d e\e% (6) 

= (/Oc*&, (7) 
{r°,P/»} = 

{p a ,P/3} = -(fi a )cB a a(3 . (8) 



In Section 2.1 we will derive the necessary and sufficient conditions for the AP bracket J5J) 
to become a Poisson bracket after an appropriate choice of quasivelocities (see Section |2F. 
For the rest of Section 2 we will concentrate on achieving that same goal by considering 
two special cases of the HPD equations ([T|-([3]): (1) where S g = {0}, known as the 
purely kinematic or nonabelian Chaplygin case, and (2) the case where Q = G, where the 
resulting equations represent a generalization of the Euler-Poincare-Suslov equations [2]. 
Let us now briefly give the details of these two special cases. 



1.1 Nonholonomic Chaplygin Systems 

Consider the subclass of nonholonomic systems with symmetry corresponding to S q = {0}, 
known as the purely kinematic case [2j [30] , where the group orbits exactly complement 
the constraints and suppose that Q ^ G. These systems are the special case of the HPD 
equations corresponding to s = (i.e. dim S q = 0, when there are no nonholonomic 
momentum equations) and are also known as nonabelian Chaplygin systems [21 [T7]. In 
the special case when Q = WL S x S r and G is either a torus action T m or acts by translations 
]R 2m , they are called abelian Chaplygin systems and correspond to the classical exposition 
of Chaplygin systems [M] where there exist local coordinates (r a ,s a ), a = l,...,n — 
2m, a = n — 2m + 1, . . . ,n such that the Lagrangian L does not depend on the s a 
coordinates, and where the constraints can be written as s a = — A^(r)r a . 

Now, from the HPD equations (|l])-([3]) we can extract the equations of motion for non- 
abelian Chaplygin systems as follows. Since s = 0, we have that l c = l(r a ,f a ,—ApfP), 
and (assuming sufficient regularity) hj^(r,p) = p a r a — l c . Then, equations ^ and 
the constraints (111) reduce to: 
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r a = (9) 

ft. = -^-0^„ do) 

T = -Al(r)r", (11) 

respectively. For easy reference later on, we also define the semi-basic two-form [13J A on 
T*M with components 

A a p(r,p) := (ti a )cB a pa , (12) 



so that the last term on the right hand side of (10) can also be expressed as A a p(dhj / [/ 'dpp) . 



Moreover, the equations of motion can be written with respect to an AP bracket: 

r a = { r a , hj^} AP , p a = {p a , hjj} AP , (13) 
where the AP bracket is the special case of ^ when ^ and ^ vanish: 

{9, k}^p ap (r,p) = {9, k} can (r,p) + A «/3^-^> ( 14 ) 

for any two functions g,k : T*M — > M, where {g, k} can (r,p) is the canonical Poisson 
bracket, 



{g, k} can (r,p) :-- 



dg dk dg dk 



dr a dp a dr a dp 



1.2 Nonholonomic Systems on Lie Groups 

Consider now another special case of the HPD equations corresponding to the setting 
where the configuration space is the Lie group G, so that there is no shape space (m = 
dim M = 0). The reduced Lagrangian becomes I = §(-/"£,£), where £ = g~ l g G Q as 
before, and / : g 1 — ^ g* is the inertia tensor. Substituting in the constraints £ b = ejQi we 
arrive at the reduced constrained Lagrangian Z c (0), and assuming sufficient regularity we 
can define the reduced constrained Hamiltonian h c (Q,p) = piVt 1 — l c . From ([!]) and Q 
the equations of motion and constraints then become: 



8h 

= -H'JcC^gf, (15) 

e = (16) 
respectively. These equations are a generalization of the Euler-Poincare-Suslov equations 

mm- 
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We can also write the equations of motion ( 15 ) as: 



Pi — {Pii h c }Ap S i (17) 

where {•, -}fp 5 is the AP bracket: 

{9,k}r P s ={Pi,Pi}^ dk 



dpi dpj 

for any two functions g,k : (g c )* — > R, where g c = {£ G = e^fi- 7 }, and where the 
bracket on the right hand side of ( [18] ) is computed by using the canonical bracket on T*G 
and then restricting to g c (see [2], Section 5.8 for more details). We note in passing that 
this bracket is merely ^ with (J7])-((8| vanishing (since m = 0). 



2 Chaplygin Hamiltonization 

To begin the generalization of Chaplygin's Theorem, we note that one can view Chap- 
lygin's time reparameterization dr = N(q)dt from the Introduction in a different way as 
follows: we have q = dq/dt = N(q)(dq/dr) =: N(q)u, which defines the quasivelocities 
w on Q (For a recent discussion of quasivelocities in nonholonomic mechanics see (6J, [8] ) . 
Thus, instead of considering which time reparameterization Hamiltonizes our system, we 
can rephrase the problem as one of finding a particular set of quasivelocities for which the 
almost-Poisson bracket ^ satisfies the Jacobi identitjj^J To that end, we need to express 
the AP bracket (tBl) in terms of the quasivelocities to, to which we now turn. 



2.1 Chaplygin Hamiltonization of the Hamilton-Poincare-d'Alembert Equa- 
tions 

Let j be the map j : (q, ou r , ujq) t— ► (q,f,Q) and define V = j*p. Then locally we have 
C c (r,u) =j*l c {r,f,ty = l c (r,r a = fu a ,Q i = fu% V a = d£ c /du r * = {dl c / dr p ){dr p / du r «) 
= fp a and similarly Vi = fpi. Moreover, define the quasivelocities uj through f = f(r, g)oo r 
and Q = f(r,g)ujQ, where / G C l is nonzero on its domain. Then we have the following 
transformation of the bracket (|5]). 

Proposition 1. Consider a nonholonomic system with symmetry governed by the HPD 
equations ^-$^>. Further, suppose that the matrix G a p := d 2 l c /dr a dr fB is invertible. 
Then the AP bracket (fBl) becomes the bracket {•, := (1/ f)j*{-, - given by: 



dG dK ~ ~ , , fdG dK dG dK 



{g, k}' m {p h v s y M d ^ d ^ + in, v a y M ^ dv ^ ^ dv 

^ ,, , dG dK dG dK\ r ~ ~ dG OK 



2 If this is successful, then some authors will also call the nonholonomic system conformally Hamiltonian (see 
Section 3.3 in [8]). 
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where 



fiWiYj* = ^k + B] 3 K (20) 

fin^hx = C^Vk + DlpK (21) 

fiPatVpYn = ElpV k + frJPv (22) 

f{r a M'M = f 5 l ( 23 ) 



and the components above are given in the Appendix by (58)- (61). 
Proof. From the reduced Lagrangian l(r,f,£), 

I = \g«pr a rP + g aa f a C + ^g^ b - V(r), 

we can form the constrained reduced Lagrangian l c (r,f,Q) by substituting in the con- 
straints (Hp in the form £ b = —A^r 13 + e b Q j . We then have p a = dl c /df a = Gapf 13 + G l a pi, 
where G af3 = g a/3 - 1g aa A% + g ab A a a A b p , G l a = (g ba - g a bA a a ) T bi and we have used e b tt j = 
T bi pi, where T ai = e"G Jt , with G 1 ^ the inverse of the matrix Gij = g a b e i e, j- Since G a p is 
invertible by assumption (denote its inverse by G a @), this leads to f 7 = G ia p a — G ia G l a Pi. 
Thus, we have: 



M^ + g^p, 



M ai {G^ a p a - G ya G i a p i ) + 9ab V b % 

M ai G^ a p a + (g ab Y u - M ai G ia G i Q ) p t 
1 



./ 



M ai G~* a V a + (g ab T u - M ai G^Gi) V % 



where we have used the definition of the V . Then, using the general relation 



(24) 



f{Pl,Vj}' M = f{f,pj}MPl ~ f{f,Pl}j4PJ + f {PhPj} 



Mi 



(25) 

which holds for all I, J — i and I, J = a, and ((6|-(|7j) along with ^ and fl24) ), we get the 
transformations in (20)-(22). Let us illustrate this for ([20|. 



From (25) we have: 



M 



f 



df_&Pj_„ _ df dpi _ 
dg a dp a % dg a dp rJ ^ j 



+ f 2 {-f{»a) c C a M e\ef) 



§^fi (ejPi - 4Pj) - ff(Va) c CZ d e b e* 

C%V h + (fC^M^G* - fK%) V k + %V a 



(26) 
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where we've used (24) in the last line of (26). Indeed, this produces (20) and (58) 



(59). The remaining transformations (21)-(22) and equations (60)-(63) follow from similar 



computations. 

Lastly, we compute j* {r ,^}^ = f{r a , Pp}' M as: 



f{r a ,VpY w = {r a Jpp} 



dr a d(fpp) dr a d(fpp) 
dr~< dp,y dp 1 dr~< 



a/3, 



which gives ( 23 ) . 



□ 



Proposition [T] gives the explicit form for the AP bracket of the HPD equations in terms of 
the quasivelocities. Now, as stated at the beginning of this section, the idea is to derive 
the conditions under which the multiplier / makes ^ into a Poisson bracket. To that 
end, we have the first main result. 

Theorem 2. Suppose that we have a nonholonomic system with symmetry satisfying 
the assumptions of Proposition [7] and let f(r,g) G C 1 be a function which is nonzero 
everywhere on its domain. Then the almost-Poisson bracket |5|) is Poisson iff f satisfies: 



B 1 - 



Am Al I Am Al _i_ Am Al 
A il A jk + A kl A ij + A jl A kv 



A kl El aB + G loP\ 



Al fik 

^ik^ja 



k/3 



fil Ak Al fik 

^ka^ij ^jk^ia 



o, 

0. 



(27) 
(28) 
(29) 
(30) 
(31) 



Moreover, the Hamiltonized equations |iJ)-Q) become, in this new Poisson bracket {-, ■} 



M ' 



r a = f{r\n M }^ V a = f{V a ,H M }^, Pi = f{Pi,HM}^, 



where TC 



M 



fh 



M- 



(32) 



Proof. In order for (19) to become a Poisson bracket it must satisfy the Jacobi identity. 



We can compute the results of this restriction, which then leads to the conditions (27)- (31 ) 



that / must satisfy. For better readability of this paper, we have left these computations 
to the second section of the Appendix. As for the second half of the Theorem, we simply 



note that since j*r = r, j*p = V and j*h 



M 



l m , then the AP bracket representation 
of the equations of motion ([l])-([3]) becomes j*(r a — {r a , hj^]j^) = 0, which gives the first 
equation in (32) since j*{-, = /{•, and similarly for the remaining. □ 



Theorem [2] represents the necessary conditions under which a given nonholonomic system 
with symmetry admitting a representation within the HPD framework can be Chaplygin 



Hamiltonized into (32). It is a generalization of Chaplygin's Theorem not only to higher 



dimensional nonholonomic systems with symmetry, but also in that it does not presup- 
pose the existence of an invariant measure. Indeed, although the conditions (27)-(31) 
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seem rather involved they are no more than a coupled set of first-order partial differen- 
tial equations in / which can be solved using any of the popular mathematical software 
packages. Thus, contrary to the traditional usage of Chaplygin's Theorem found in the 
literature (where the reducing multiplier is typically guessed at by knowing the system's 
invariant measure and degrees of freedom), Theorem [2] does not require the invariant 
measure and eliminates the guesswork. In fact, as we will show below, for nonabelian 
Chaplygin systems we will recover the second part of Chaplygin's Theorem from the 
analogous / conditions in that case, showing that for these types of nonholonomic sys- 
tems it is more advantageous to solve the corresponding / conditions instead of guessing 
at the multiplier using the invariant measure density, since if a solution exists then we 
will get the invariant measure density for free. Moreover, as we will also show below 
and have already mentioned above, one can now also apply Chaplygin Hamiltonization to 
systems for which both Chaplygin's Theorem and the guesswork above are inapplicable: 
nonholonomic systems which do not possess an invariant measure density. The Chaplygin 



sleigh (which we discuss in Section 4.5) is perhaps the best illustration. 

Let us now turn to the Chaplygin Hamiltonization of the two special cases considered in 
Sections 1.1 and 1.2 above. Indeed, consider now the special cases of the HPD equations 
given by the nonabelian Chaplygin nonholonomic systems, as in Section 1.1, and the 
Euler-Poincare-Suslov case of Section 1.2, where there is no shape space. From Theorem 
[TJ in these cases the quasivelocity transformations read r a = f(r)u r a (we shall write 
this simply as r = f(r)u henceforth) and Q l = f{g)ujQ i i (we shall write this simply as 
Q = f(g)u henceforth), respectively, where / G C l is nonzero on its domain. Moreover, 
defining the maps j r : (q, u) i— ► (q, r) and jn : (u) i— > (Q) we have the following Corollary 
of Theorem [2] 

Corollary 3.(1) For a nonabelian Chaplygin nonholonomic system (L, G, T>) described by 



Q)-( li) with almost-Pois son formulation (13)-(14) satisfying the assumptions of Propo- 
sition^ the necessary and sufficient conditions for Chaplygin Hamiltonization (using 
dr = fdt ) on j*M are that f satisfy 



f r {9,k} c A h P ap (r,p) = f{G } K} can (r,V), 



(33) 



for all a, v, 8 = 1, . . . , m (recall from Section 1 that m = n — k = dim(M) ), and where 
{•) - }can is the canonical bracket on j*M. 



(2) For an Euler-Poincare-Suslov nonholonomic system described by (15)- (16) with almost- 
Poisson formulation (11)-(18) satisfying the assumptions of Proposition [i| the necessary 
and sufficient conditions for Chaplygin Hamiltonization (using dr = fdt) on Jq(0 c )* are 
that f satisfy 



Jn{g,k} E A P P S = f{G,Ky} (0C) \ 
where {•, •}'!" < ' 9 ^ is the (minus) Lie-Poisson bracket on j^(g c )*. 



(34) 



Equivalently, in local form conditions (33) and (34) read 
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9f_ r df n 9f 

Q r 5 ^ al/ Q r u ^ aS Z Q r a USl - 



^km^ij + ^im^ki + ^im^ik ~ ^> h I 



} jrn Lj ki 



m 
im u jk 



l,...,n 



(35) 
(36) 



respectively, where S l km 



km 



Proof. (1) Let us consider the nonabelian Chaplygin case first. Since this case corresponds 
to the situation in which s = dim S q = 0, then only Greek indices survive in Theorem 
[2j Moreover, using that in this case g a b = g aa = for all a, a and / is independent 
of g, we can extract the relevant Hamiltonization conditions from (27)-(31). The only 



non- vacuous condition amongst those in Theorem |2| is then = 0. From (|63|) this leads 



to the condition (33), and its local form in (35). 



(2) For the Euler-Poincare-Suslov case, since this corresponds to the special case of the 
HPD equations in which m = dim M = 0, then only the Latin indices survive in Theorem 
[2j Thus, g a p = g aa = for all a, a, (3, and since / is independent of r, the only non-vacuous 
condition amongst those in Theorem [2] is condition (28). However, note that because of 
the fact that m = 0, the first term in ( J58| ) vanishes, giving the condition (34), and its 
local form in (36). □ 



For completeness, we should note that based on the results of Corollary [3] we can write 
the reduced constrained mechanics of a nonabelian Chaplygin and Euler-Poincare-Suslov 



nonholonomic system from (32) as 



fi^^Mjcan and Vp = f{V/3,Hj4} c 



(37) 
(38) 



respectively. The reader familiar with the usual treatment of Chaplygin's theorem will 



note the absence of the reparameterized time r in (37). In fact, in the context of Chap 



lygin's work, as well as to compare directly with [7j, we note that the quasi-Hamiltonian 
forms (37)-(38) are the "t-time" analogues of the Hamiltonian forms stated in the classical 
Chaplygin Reducibility Theorem in "r-time," and the two are related through r = fr', 

V = fV', where r' = dr/dr and V' = dV/dr. 



Now, given the more general conditions in Corollary [3] (which are valid for nonabelian 
Chaplygin and Euler-Poincare-Suslov systems in arbitrary degrees of freedom), let us 
proceed to extract both parts of Chaplygin's Reducing Multiplier Theorem as special 



cases. 



2.2 Chaplygin's Reducing Multiplier Theorem 

We now specialize to the case when m = dim M = 2 (the two degree of freedom case) to 
extract the first part of Chaplygin's Theorem. 
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Corollary 4. The necessary and sufficient condition for a Chaplygin nonholonomic sys- 
tem (L,G,T>) in two degrees of freedom (m = 2) to be Chaplygin Hamiltonizable is that 



dKL OK 2 . . 

L - (39) 



equivalently that the system |g|)-(|70|) possess a nonzero invariant measure density 
N(r) G C 1 . The multiplier is then given by f(r) = ef K i2 dr2 = N. 



Proof. From Corollary |3| the only independent conditions in (33) in the two degree of 
freedom case are: 



df_ 
dr 2 



@-/jfi)a«-(^+/J4)G- = o. 

Since we have assumed that G a p is invertible, the necessary and sufficient condition for 
the satisfaction of these equations is that the parenthetical terms vanish. The resulting 



set of equations is soluble iff (39) is satisfied, in which case / is given in explicit form as 



in the Corollary. However, as we showed in [23], if the the constrained reduced system 
(gl-gOj has an invariant measure, then its density N (for m = 2) satisfies: 



K{ 2 



1 dN 

Ndr 2 ' 



K 2 

-"■12 



1 dN 
Ndr 1 ' 



One sees immediately that this satisfies (39), and hence / = N is a multiplier 



□ 



Moving on to the second part of Chaplygin's Theorem, the Proposition below shows that 



it too follows from the Hamiltonization condition (33). 



Proposition 5. Suppose f satisfies the conditions of Corollary [3| Then the original 
system (^)-(lO) has an invariant measure with density f Tl 



Cm—l 



Proof. Suppose / satisfies (35). Multiplying by r and f v and adding results in: 



1 

7 



df 



QyP^ a 3r a 
Comparing the coefficients yields: 



pp r p 



A 



Pa 



df_ 

f \dr? 



Pc 



df 
dr' 



-PP 



(40) 



where we remind the reader of the definition of A from (12), and have used p Q = dl c /df 



(we will drop the tildes in p here). We also note that the relationship (40) was also 
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presented as a sufficient condition for the existence of an invariant measure by [36J (see 
also [17J), but here is derived from the conditions of Corollary |3} Thus, for a Chaplygin 
Hamiltonizable system the second term on the right hand side of (10) can be written in 
terms of / as in (40). 

Now, suppose X n h = r a d r a +p a d p a is the nonholonomic vector field solution to the system 
@-(|ll}. We will show that f m ~ l is an invariant measure density by showing that the 
vector field f m ~ 1 X n h has zero divergence. A straightforward calculation yields 



div U m ~ l Xnh) 



dr c 



dp c 



flu -2 ■■xit 1\ ®f i "> 



(41) 



and a simple calculation of the last term in (41) using (40) then shows that the divergence 
does indeed vanish and completes the proof. □ 

Corollary [3] yields the necessary conditions for the Chaplygin Hamiltonization of the 
nonabelian Chaplygin system ([£)-( 11), which locally are the first-order partial differential 
equations (35) in r. Proposition 5] then completes the generalization by providing us with 



the invariant measure density given a solution to (33). We should stress, however, that 
the converse of Proposition [5] is not true in general. That is, given an m > 2 degree of 
freedom nonholonomic system with an invariant measure, its Chaplygin multiplier / may 
or may not coincide with the invariant measure density (or any other smooth function of 
it). This is most easily seen by first assuming that the m degree of freedom Chaplygin 
system has an invariant measure with density / m_1 , so that the right hand side of (41) 



vanishes. Inserting the resulting equation for df/dr a into (35) yields 



2G S „K. 



v3 



[171 



1) (K S G,u 



Kt.G 



fj,8) 



(42) 



The conditions in (42) are, for m > 2, conditions solely arising from the nonholonomic 
system itself, as (|42[) depends only on the metric of the Lagrangian and the curvature 
of the connection (for m = 2 (42) is vacuous, a manifestation of Corollary BR Thus, if 



the metric and curvature of the connection of a Chaplygin nonholonomic system interact 



precisely as in (42) then the converse of Proposition 5 holds. However, given the rarity of 



such an event, we believe that the ordering of the Hamiltonization process for nonabelian 



Chaplygin systems that first begins by attempting to solve the conditions (35) and then 
extracting the invariant measure density from Proposition [5] is bestQ 

In another direction, it may be the case, however, that (|35l) does not have a solution. This 



does not mean that the system is not Chaplygin Hamiltonizable though, since it may still 
possess more symmetries which further reduce the degrees of freedom, and which allow 
one to seek such a solution on the second reduced phase space. We illustrate such a 
situation in the next section, making use of the classical Routhian [5] to explore the effect 
of additional simple symmetries in (|9])-(10) on its Chaplygin Hamiltonizability. 



3 Moreover, it should also be clear that for m > 2, using any smooth function F(f; m) will also lead to conditions 
similar to ( 42 l . 



4 Nonetheless, it is impressive that some authors |19l 1201 [21] have effectively found nonholonomic systems for 



which (421 is satisfied. 
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2.3 Momentum Conservation and Chaplygin Hamiltonization 



Suppose that (35) has no solutions, but that the nonholonomic system possesses momen- 
tum conservation laws that we have yet to account for. With the aid of these conservation 
laws, we can apply the reduction process to further reduce the degrees of freedom of the 
system and re-attempt a Hamiltonization on the second reduced space. 

In order to illustrate this in a simple manner, we restrict ourselves in this section to 
Chaplygin systems which we will call nonholonomic cylic. By this we mean that we have 
an abelian Lie group H acting on M = Q/G by M 3 r a = (w a ,v l ) t— > (w a ,v l + k 1 ), 
h G H, where i — 1, . . . I — dim(H) and a' = 1, . . . , m — I and such that A a >i = Vi, a' 
(we shall hereafter denote the nonconserved conjugate variable indices w with a prime) 
and such that the action leaves the Lagrangian and constraints invariant. Under these 
assumptions the v l equations in ( 10 ) lead to the momentum conservation law^jand we can 



thus set pi = Aj = constant and perform a partial Legendre transform in the v 1 variables 
to form, analogous to the classical Routhian [321 [33] , the constrained Routhian R c (w,w) 
defined by 

R c (w,w) := [l c (w,w,v) - XiV l ] p . =Xi . (43) 

Now, using the well-known fact [32] that the Euler-Lagrange expressions of the noncon- 
served variables of l c are equivalent to the Euler-Lagrange expressions of the nonconserved 
variables of R c , we can write the Lagrange-d'Alembert equations for R c as 

d L dR^_dR^__(dl \ ™ .p> (u) 

along with the conservation equations pi = 0. Furthermore, the last term on the right 



hand side of (44) can be rewritten in terms of the Routhian: 



31 

lvlaa(j a ' 13 ' 



31 

M aa G^-^ + M aa G a % 
3w e 



3R 

, 3R r 



pa 



(45) 
(46) 



where G'' 



G m — G ae 'G^ 3 Gj e i, assuming the invertibility of C? 



(3 2 l c /dv l 3v J ) as well 



as that of the kinetic energy matrix of l c (we also remind the reader of the definition of 
the K^q in (64)). Moreover, we shall henceforth denote the parenthetical term in (45) by 

(¥R c y. 

5 Since the fMnvariance implies that l c does not depend explicitly on the v\ these variables are cyclic and 
produce momentum conservation laws in unconstrained systems. However, due to the presence of the A a /3, 
cyclic variables are not enough to produce the conservation laws, hence the introduction of the terminology 
"nonholonomic cyclic." 
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Then, since R c can now be interpreted as a function on T(M/H), we can now attempt 
to Hamiltonize (44) on this second reduced space. To that end, defining the maps j : 



(w,w) and j : (r, Aj) i— > (r,pi), we have the second main result. 



Theorem 6. Suppose that the nonabelian Chaplygin nonholonomic system given by Q)- 
(11) is not Hamiltonizable by Corollary^ but is nonholonomic cyclic. Further, suppose 

that the kinetic energy matrix ofl c and the sub-matrix ~G l ° := [d 2 l c /dv l dv^) are invertible. 
Then if there exists a multiplier f{w), nonzero everywhere on its domain with f(w) G C 1 , 
satisfying 



dK c 
duj 



((WK c )',B(uJd w )), 



(47) 



where IZJw, u 



j R c and (F7Z C )' = j (¥R C )' , the reduced system is Chaplygin Hamil- 
tonizable on M' = M/H under the choice of quasivelocity w = fu. Furthermore, assum- 
ing sufficient regularity its dynamics on T(M/H) can be written in the quasi- Hamiltonian 
form 



w ' 



f{w a ,Hm'}ap, Vp> = f{P > 8/,'HM'}AP, 



(4f 



where V' a , = dlZ c /duj a ' and Hm' = ^'^'a' ~ H'clu-tp' ^ s ^ie Hamiltonian and the almost- 
Poisson bracket is defined by: 



{G,K} AP {w',V') = {G, K} can (w', V') — fK l a , ,\, 



dG dK 



(49) 



where {-,-} C an is the canonical bracket on T(M/H). Moreover, the bracket automatically 
satisfies the Jacobi identity for dim(M') = 2. 



Proof. Under the quasivelocity transformation w = f(w)u the constrained reduced equa- 



tions (44) become (taking into account (46)) 



d dTZ c dTZ c 

Zr J ~ 



dt du c 



dw a ' 



, dTZ c nA 2 * 



(50) 



where W*, 3 , := fK%, a , — C^,g,. Now, if / is chosen to satisfy (47), then the parenthetical 
term in (50) vanishes. By defining the Hamiltonian 7~Cm> as m the statement of the 



Theorem, the equations of motion can then be written as in (48) with the almost-Poisson 



bracket (|49j). Lastly, a straightforward computation shows that the Jacobi identity is 
automatically satisfied for dim(M') = 2, owing to the fact that the non-canonical part in 
(49), {Va'i V/3'}ap, is independent of the momenta. □ 



Theorem [6] will be used below when discussing the Chaplygin sphere (a classic example of 
how the apparent failure of Hamiltonizability can be reversed in the presence of momentum 
conservation laws) and the Snakeboard as well. However, let us remark here that in [7] 
the authors consider an extension of the Chaplygin method to the case where gyroscopic 
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forces are involved, analyzing the Hamiltonization of nonholonomic systems in two degrees 
of freedom. In their subsequent bracket description of the Hamiltonized mechanics there 
appear non-canonical parts which in this paper manifest themselves as the second term 



in (49), resulting from the second reduction to M/H . Indeed, we can compare the results 



of Theorem ^ to the exposition in [7J by first noting that the non-canonical part of the 



bracket (49) is the many degree of freedom analogue of S in equation (3) of [7j. This is 
best seen by defining S a >pi := —fK l a ,a,\i along with the 2- form Q := S a >[3>dw a A dw& . 
The two-form Q is exact when dim M = 2, and in that case (or any other case when it is 
exact), one can then locally write Q = d(3, where (3 = W a i{w)dw a , and as the authors in 



[7] point out, the constrained reduced equations (48) can then be rewritten as 



d dK w _ dTl w 

dtdcu*' 1 dw<*' ' 1 J 

where Rw{w,uj) = R c {w,uj) + W a >(w)uj a ' . 

As a preliminary application of the results above, we shall now use Corollary [3] to extend 
prior Hamiltonization results from [23J of a class of nonholonomic systems known as con- 
ditionally variational systems, which are nonholonomic systems that can be Hamiltonized 
in full subject to the imposition of initial conditions that satisfy the constraints. 



3 Conditionally Variational Systems in the Quasivelocity Con- 
text 

In [23] we discussed the notion of a conditionally variational nonholonomic system. 
Briefly, these systems have the property that the constrained Euler-Lagrange equations are 
Lagrangian, hence making it possible to express the constrained dynamics in a variational 
manner. However, we showed that under certain additional requirements for the original 
system's Lagrangian one can construct the "variational" Lagrangian Ly whose Euler- 
Lagrange equations reproduce the nonholonomic equations when the initial conditions 
are chosen to satisfy the constraints (as they must anyway). Hence, such nonholonomic 
systems can be realized as variational systems provided the initial conditions satisfy the 
constraints and that the original Lagrangian satisfies certain requirements. 

In that paper we showed that such a realization was possible only in the cases when the 
nonholonomic system possessed an invariant measure with constant density N{r) (the 



vertical rolling disk of Section 4.1 is such a system). However, using Corollary |3j we 
can now extend the results in |23j to a more general setting if we instead focus on the 
Chaplygin Hamiltonized system. To that end we have the following result: 

Theorem 7. Suppose that for a given abelian Chaplygin nonholonomic system (L, G, V) 
with constraints given by 

<P a (q,q) = s a + A a a (r)r a , a = 1, . . . , k < n, 

where q = (r, s), we have found an f as in Corollary^ above and let £(q,tu) := L(q,f = 
fu, s = fuj) and <fi a (q,u) = (f> a (q,r = fu, s = foj). Then if the matrixg a b := {d 2 C/ dui a dui h ) 
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is invertible, the nonholonomic mechanics of the original system can be derived from the 
(almost) Euler-Lagrange equations 



d dC v dC v 

f 7^T = ' I=l,--,n, (52) 



dt duo 1 dq 
by using the Lagrangian £y{q,uj) defined by 

1 dC 

C v (q, u) = C(q, uj) - -J-Q^^ili w )> (53) 
and imposing the nonholonomic constraints initially. 

Proof. The existence of an / which Hamiltonizes the system guarantees, by part (2) of 
Proposition 3 in [23], that the system (£(q,ui),(f>(q,ui)) is conditionally variational after 
the reparameterization dr = f(r)dt. Then, the Theorem follows by Proposition 5 of [23] 
again. □ 

Theorem [7] extends Chaplygin's Theorem in a different direction. Unlike Corollary |3j it 
gives one a method to Hamiltonize the entire system (similar to some of our earlier work 
[I]) after a time reparameterization. Although it is impossible for a nonholonomic system 
to be Hamiltonian [2] , Theorem [7] begins to answer the open problem briefly discussed 
in [8] and elsewhere of lifting the Hamiltonization of the reduced problem to the whole 
system. We illustrate this and our other results below by applying these ideas to some 
common and well-known nonholonomic systems. 



4 Examples 

The simplest illustrations of the above results can be found in low dimensions, specifically 
the two degree of freedom case. Although this is the original setting for Chaplygin's 
Reducibility Theorem, we will discuss the Chaplygin sleigh (which, due to its lack of an 
invariant measure, cannot be handled by Chaplygin's Theorem), among other things, and 
also take this opportunity to illustrate Theorem [7] as well. 



4.1 The Vertical Rolling Disk 

Consider the nonholonomic vertical rolling disk pictured in Figure 1 below with configu- 
ration space Q = M. 2 x S 1 x S 1 and parameterized by the coordinates (x,y,9,(p), where 
(x, y) is the position of the center of mass of the disk, 9 is the angle that a point fixed on 
the disk makes with respect to the vertical, and <p is measured from the positive x-axis. 
This system has Lagrangian and constraints given by: 

L = X -m{x 2 + f) + \l0 2 + \ J<p\ 
(f) 1 = x — Rcos tpO = 0, 

<p 2 = y- Rsmcp6 = 0, (54) 
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where m is the mass of the disk, R is its radius, and /, J are the moments of inertia about 
the axis perpendicular to the plane of the disk, and about the axis in the plane of the 
disk, respectively. 




Figure 1: The Vertically Rolling Disk. 



This system has an invariant measure with constant density N, which we can take without 
loss of generality to be unity. Corollary [4] applies and since f = N = const, then we can 
simply take q = u so that the new quasivelocities are merely the original g's. A short 
computation then shows that the system (54) satisfies Theorem [7j and the variational 
Lagrangian is computed through (53) to be 



Lv(q,q) 



-m(x 2 + y 2 ) + -AQ 2 + ~^J^ 2 + mR9(xcos ip + ysin </?). (55) 



The Lagrangian (55) first appeared in [23], and a short computation shows that apply- 
ing the initial conditions 1 (O) = 0, <f) 2 (0) = to the Euler-Lagrange equations for Ly 



reproduces the nonholonomic equations for the system (54). 



This simple example illustrates the case when Hamiltonization is automatic (i.e. / 
const) and thus the system is conditionally variational as well. 



4.2 The Nonholonomic Free Particle 

Consider a nonholonomically constrained free particle with unit mass (more details can 
be found in [2]), and Lagrangian and constraint given by 



L=l(x 2 + y 2 + z 2 ) } 

0(?> q) = z + xy = 0. (56) 

The system possesses an invariant measure with density N(x) = (1 + x 2 ) -1 / 2 , and thus 
by Corollary [4] the system is Chaplygin Hamiltonizable with f(x) = (1 + x 2 ) -1 ^ 2 and 
quasivelocities defined by uj = y/l + x 2 r, where r = (x, y). 

To illustrate Theorem [71 note that C(q,u) = (1/2)/ 2 (c<j 2 + lo 2 + uj 2 ) and that since g zz = 
(1/2) f 2 , Theorem u\ applies and we have Ly given by: 
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2(1 + x 2 



2XUJ X UJyj 



Then, computing equations (52) gives: 



Co,, 



'l + x 2 fl 2 ' 



d 
db 



{f 2 {uJ Z + XUJy)) 



0. 
0. 



(57) 



Now, the last line of (57) reads (d/dt)(f(f)(q,uj)) 

f(q(t))<t>(q(t)Mt)) = /(g(0))#7(0)M0)). 
then <fi(q(0),u(0)) = 0, and hence <f)(q(t),u)(t)) 



0, which gives the conservation law 
Thus, if the constraints are satisfied initially, 
for all t. Recalling that q = fu, this 
then expresses the conservation in time of the original constraint equation (56). After 



imposing the constraints initially, one can then use the quasivelocity definitions to then 
transform Co — > r and recover the original nonholonomic mechanics that results from the 



application of the Lagrange-d'Alembert principle to the system (56). Thus, although (52) 
is not Hamiltonian, as has been the theme in this paper, it is after Chaplygin's time 
reparameterization (and the imposition of initial conditions satisfying the constraints). 
Thus the nonholonomic free particle, like the vertical disk, is Hamiltonizable but since 
/ 7^ const it is only conditionally variational after a reparameterization of time. 



4.3 The Chaplygin Sphere 

The Chaplygin sphere is a sphere rolling without slipping on a horizontal plane (see [2J [7]) 
whose center of mass is at the geometric center, but the principal moments of inertia are 
distinct. In Euler angles (9,ip,ip) the Lagrangian and constraints are: 



— # sin (p + tp cos (p sin 6j 
^(V + </>cos#) +-(x 2 + y 2 ), 

(j) 1 = x — 6 sin tp + ip cos ip sin 9 = 0, 
<p 2 = y + 9 cos ip + if sin ip sin 9 = 0. 

where Jj are the moments of inertia about the center and where we have assumed the ball 
to have unit radius and mass. 

Since q = (x,y,9,ip,(p) and the constraints and Lagrangian are cyclic in x,y, we can 
consider this to be an abelian Chaplygin system. The system has an invariant measure 
whose density N(9,<p) is in general non-constant [7]. 



h (• • \ 2 

— [9 cos (p + ip sin ip sin 9j + 



+ 
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Applying Corollary [3] shows that there does not exist an / which Hamiltonizes the three 
degree of freedom base dynamics given by (|9j)-(10) when viewed as an abelian Chaplygin 



system. However, it is easily seen that ip is a nonholonomic cyclic variable and leads to 
the momentum conservation law p^ = Xf. Thus we can form the constrained Routhian 
as in (43) and further reduce the dynamics to M' = S 1 x S 1 . We can then Hamiltonize 
on M' through Theorem|6l from which (47) shows that / = N(0,ip). The non-canonical 



part of the almost-Poisson bracket (49) is then computed to be 



{V'i,V' 2 } = -\4h + l)f sme{hcos 2 if + I 2 sm 2 ip + 1), 

and by the same Theorem since dim M' = 2 we know that this bracket satisfies the Jacobi 
identity and hence is indeed a Poisson bracket. This matches the result obtained in [7j 
and is an example of a system that although is not Hamiltonizable at when viewed as a 
three degree of freedom abelian Chaplygin system is in fact Hamiltonizable on the second 
reduced space M' of dimension 2. Moreover, it also serves to illustrate the discussion at 
the end of Section 12.31 



4.4 The Snakeboard 

Another example of Theorem [6j whose greater importance we will discuss in the Conclu- 
sion, is the Snakeboard [2J [31] . This system is modeled as a rigid body (the board) with 
two sets of independent actuated wheels, one on each end of the board. The human rider 
is modeled as a momentum wheel which sits in the middle of the board and is allowed to 
spin about the vertical axis, see Figure 2. 




Figure 2: The Snakeboard. 

The configuration space is Q = SE(2) x S 1 x S 1 and the Lagrangian L : TQ — > R and 
constraints are given by: 



x z + y z + r + tp z + 2ipe + 

x = — cot cos 99, 
y = — cot cf) sin 99, 
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where we have set the mass m, moments of inertia, and the distance r from the center of 
the board to its wheels equal to unity. Here (x, y, 9) represent the position and orientation 
of the center of the board, ip the angle of the momentum wheel relative to the board and 
0i and 02 the angles of the back and front wheels relative to the board. Here we've made 
the simplification that 0i = —02, as in [21 [3Tj . 

As stated, we can view this system as an abelian Chaplygin nonholonomic system with 
three degrees of freedom. Its equations of motion are given in [3T] as: 

Pe = ~ - sec esc 4>(p e - p^)p^ 9 
]ty = 0, 
Pi> = 0, ip 



tan 2 (j)(p e - p^) 
1 



2 ^ 

P^ 



sin 



cos^ 



Since this system satisfies the conditions of Theorem [6] we can set p^ — = const. 
and focus on Hamiltonizing the reduced system. Given that this reduced system has the 
invariant measure iV(0) = tan0, which is independent of tp, by Corollary 111 / = N. The 



non-canonical part of the almost-Poisson bracket (49) is then computed to be: 



{V[ } V' 2 } = sec 2 ( f ) X lP , 

and by the same Theorem we know that this bracket satisfies the Jacobi identity (since 
the reduced system has two degrees of freedom) and is thus a Poisson bracket. 



4.5 The Chaplygin Sleigh 



The Chaplygin Sleigh [21 El [HI HH E3 EH] consists of a rigid body in the plane which is 
supported at three points, two of which slide freely without friction while the third is a 
knife edge, a constraint that allows no motion perpendicular to its edge. The configuration 
manifold Q = IR 2 x S 1 , where (x, y) are the coordinates of the contact point while 9 is the 
angle the knife edge makes with the x-axis, see Figure 3 below. Moreover, we suppose 
here that the center of mass of the system C is not on top of the knife edge (if it is, then 
one can show [2] that the sleigh reduces to another nonholonomic system known as the 
knife edge, which possesses an invariant measure). 

The Lagrangian L and constraints are given by: 



L = - (i 2 + y 2 + 29 2 - 2(x sin^ + y cos 9)9 j , 
y cos 9 — x sin 9 = 0, 

where for simplicity we have set all parameters to unity. Since the Lagrangian and con- 
straint are left invariant on the Lie group G = SE(2) we can treat the problem within the 
Euler-Poincare-Suslov framework. Defining £ = g~ l g, where g = (x,y,9), we can write 
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Figure 3: The Chaplygin Sleigh. 



the Lagrangian L in terms of £ as /(£) = £3 + (1/2) (£1 + £| ) + £2^3, and the constraint as 

e = 0. 

With the structure constants given by Cf 3 = —1 = —C\ z and all other zero we see that 
/ = const, satisfies Corollary [3j which agrees with the recent result of |25j . 

This system is of critical importance in the study of Hamiltonization since unlike Propo- 
sition [5j the Chaplygin sleigh shows that just because a system is Hamiltonizable does 
not imply that it possesses an invariant measure. Indeed, although it is well-known the 
Chaplygin sleigh does not possess an invariant measure (21 [25], as we've seen above this 
system is nonetheless Hamiltonizable. Thus, unlike for the nonabelian Chaplygin case, the 
Hamiltonizability of Euler-Poincare-Suslov systems does not automatically imply that the 
system possesses an invariant measure, and thus Chaplygin's Reducibility Theorem be- 
comes inapplicable (due to the non-existence of an invariant measure)^} However, thanks 
to the results of Corollary |3j we may still be able to Hamiltonize, or, more properly, 
"Poissonize" (which would be the better term here since there isn't an invariant measure, 
as discussed in the Introduction). 

4.6 A Mathematical Example 

Consider the following mathematical example due to Iliyev [27] . The Lagrangian and 
constraints are given by: 

q 4 = g 2 tan(g x ), 
q 5 = g 3 tan(g x ). 

This is a nonholonomic system with three degrees of freedom (m — 3), and thus Chaply- 
gin's Theorem is inapplicable. Within our framework, it can most easily be treated as an 

6 However, it is interesting to note that Chaplygin did apply his reducing multiplier method to the Chaplygin 
sleigh in [15], but only after introducing "quasicoordinates." 
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abelian Chaplygin system. Solving the conditions in (35) in MAPLE yields / = cos(g 1 ). 
Moreover, as a check of Proposition |5j one can show that the system's invariant measure 
density is N = cos 2 (g 1 ), which indeed is equal to / m_1 , as the Proposition suggests. 



5 Conclusion and Future Directions 

Chaplygin's Reducing Multiplier Theorem has long allowed an interesting investigation 
of some nonholonomic systems in terms of the quasi-periodic orbits that result from 
the consideration of the time reparameterization dr = fdt and the Hamiltonian-like 
structure it produces. Perhaps because of this, and its success in studying nonholonomic 
systems using methods from unconstrained mechanics, it has attracted much attention 
in the recent decades as interest in nonholonomic systems has grown. However, as we 
have mentioned, the use of Chaplygin's results, and much of the subsequent research 
that has followed it, has been confined to systems possessing an invariant measure and, 
typically, also in two degrees of freedom. In addition, and partly due to these confines, 
Hamiltonizability of nonholonomic systems with symmetry in arbitrary dimensions has 
remained untouched (with the only results [21] and [28] known to the authors arising by 
construction) . 

The present work addresses these two main aspects of Chaplygin's work, extending the 
results to nonholonomic systems in arbitrary degrees of freedom not necessarily possessing 
an invariant measure (a central assumption of the research on Chaplygin's work to date). 
As such, this latter result alone represents a possibly new direction for the study of the 
integrability of nonholonomic systems, with a Hamilton- Jacobi theory based on it now 
possible (it would be interesting to develop this and compare it to the Hamilton- Jacobi 
theory of nonholonomic systems recently presented in [26j. In fact, as we have mentioned, 
Chaplygin's Theorem was used in conjunction with the Hamilton- Jacobi method implicitly 



in [151 EI], albeit in "quasicoordinates."). The local conditions (27)-(31), or their special 



cases (35) and (36), also enable the search for Hamiltonizable nonholonomic systems to be 
converted into the search for solutions to certain partial differential equations, a task which 
can be considerably simplified by making use of any of today's mathematical software 
packages and which eliminates the guesswork involved in current uses of Chaplygin's 
Theorem. 

We should also point out the interesting role that symmetry plays in Chaplygin Hamil- 
tonization. For example, considering the Chaplygin sleigh as an abelian Chaplygin system 



it is immediately seen to be impossible to Chaplygin Hamiltonize (the equations (35) have 
no solution), yet as we showed in Section 4.5 it is Chaplygin Hamiltonizable when consid- 
ered as an Euler-Poincare-Suslov system. This suggests to us that the choice of symmetry 
group affects the Hamiltonizability of the system in question. We expect to pursue these 
issues in future research. 

Finally, we note that the multi-dimensional Veselova system and multi-dimensional Chap- 
lygin sphere have recently been Hamiltonized in [21] and [28], respectively. However, 
the methods and conditions for Hamiltonization presented here are inapplicable to those 
Hamiltonizations due to the particular Hamiltonization methods used by the authors. 
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In the former, the authors constructed redundant coordinates and showed that the so- 
lutions of the multi-dimensional Veselova system can be mapped isomorphically into the 
solutions of an associated different Hamiltonian system known as the Neumann system. 
Within the framework of the methods presented here and in our previous research [3] , this 
would be equivalent to the statement that after an appropriate time reparameterization, 
applying the inverse problem of the calculus of variations to the resulting system would 
yield the Neumann Lagrangian as a solution. 

In the latter case, the author Hamiltonizes the multi-dimensional Chaplygin sphere by 
constructing redundant coordinates and effecting a time-reparameterization. It is then 
shown that the reduced mechanics of the higher dimensional nonholonomic Chaplygin 
sphere emerge as the restriction to the invariant submanifolds of the Hamiltonian system 
resulting from the time reparameterization. In our previous research [1] we called these 
type of systems associated second-order systems. However, the main difference between 
our work there and the construction in [28J is that we constructed associated second-order 
systems for the original nonholonomic system (not the time reparameterized one). 

Given the above discussion, we therefore expect that the aforementioned multi-dimensional 
Hamiltonizations can be realized as special cases of a synthesis of the general (yet mostly 
disjoint) methods presented here and in earlier work jlj (see also [22J). 
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7 Appendix 

7.1 The Components of the Quasivelocity AP Bracket 



The components of (20)-(22) are given by: 
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A 



C 



k 

ia 



T?k 

r a/3 



fCt d e\e]M a ,G^G k a - [fK% - C* 

-fC a cd e^M a ,G^, 
df df a 



-g° d A d a + fF? a (g ab T bk - M a ,G^G k ) 



At = fF? a M a ,G^-% 9d -ei8i 



dg 

fB b ap (M h G^G k e - g bd T dk ) , 

{A h J} - A b ,5l) + (fKj a - Cjg) 



(58) 
(59) 

(60) 

(61) 
(62) 
(63) 



where 



with M„ 



g ab A b c 



21 /3a 

K k 

j* 



C 



M bt G el Bp aJ 
n l r a p c P d V bk 

gab { ^cd e i e j L i 
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df 
dg 



~ a 9l m ~ efS k ) 



(64) 
(65) 

(66) 
(67) 



7.2 Calculation of the Jacobi Identity for the Quasivelocity AP Bracket 



Here we illustrate the calculation of the Jacobi identity for the bracket (19). Since it 



is well-known [2] that the Jacobi identity is satisfied iff it is satisfied for the component 
functions, we need only calculate it for all combinations of x = (r 7 , Vi, V a ), i.e. we require 
{x 1 , {x J , xK }'jy}jj + cyclic = for all I,J,K= (a, i, a). 

As an example, consider x = (r 7 , V a , Vp). Then we have: 



{r\{V a ,VpY IT y w + cyli C = 0, 

0, 



J_A7 P 



which gives the third equation in (27). Similarly, considering x = (r 7 ',Vi,Vp) gives: 



{r\{v l ,Vpy w y w + cj\i C = o, 
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which gives the second equation in (27). Similar computations lead to the remaining 
conditions in (27). 
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